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ABSTRACT 

Context. The interaction of microquasar jets with their environment can produce non-thermal radiation as is the case for extragalactic 
outflows impacting on their surroundings. Observational evidences of jet/medium interactions in galactic microquasars have been 
collected in the past years, although few theoretical work has been done regarding the resulting non-thermal emission. 
Aims. In this work we investigate the non-thermal emission produced in the interaction between microquasar jets and their environ- 
ment, and the physical conditions for its production. 

Methods. We have developed an analytical model based on those successfully applied to extragalactic sources. The jet is taken to be a 
supersonic and mildly relativistic hydrodynamical outflow. We focus on the jet/shocked medium structure when being in its adiabatic 
phase, and assume that it grows in a self-similar way. We calculate the fluxes and spectra of the radiation produced via synchrotron, 
Inverse Compton and relativistic Bremsstrahlung processes by electrons accelerated in strong shocks. A hydrodynamical simulation is 
also performed to further investigate the jet interaction with the environment and check the physical parameters used in the analytical 
model. 

Results. For reasonable values of the magnetic field, and using typical values for the external matter density, the non-thermal particles 
could produce significant amounts of radiation at different wavelengths, although they do not cool mainly via radiative channels but 
through adiabatic losses. The physical conditions of the analytical jet/medium interaction model are consistent with those found in 
the hydrodynamical simulation. 

Conclusions. Microquasar jet termination regions could be detectable at radio wavelengths for current instruments sensitive to ~ ar- 
cminute scales. At X-rays, the expected luminosities are moderate, although the emitter is more compact than the radio one. The 
source may be detectable by XMM-Newton or Chandra, with 1-10 arcseconds of angular resolution. The radiation at gamma-ray 
energies may be within the detection limits of the next generation of satellite and ground-based instruments. 

Key words. ISM: jets and outflows - X-rays: binaries - Radiation mechanisms: non-thermal 



1. Introduction 

Microquasars (MQ) are radio emitting X-ray binaries (REXB) 
that display relativistic ejections (Mirabel & Rodriguez 1999). 
The binary system is formed by a compact object, a neutron 
star or a black hole, and a normal (non-degenerate) star. The 
relativistic ejections may be transient or persistent (Ribo 2005 ), 
with duty cycles of jet activity that change from source to source. 
The steady jets are formed during the so-called Low/Hard state 
(Fender et al. 2001), and are expected to be only mildly rela- 
tivistic (Gallo et al. 2003). Presently, about 15 MQs are known 
in the Milky Way (Paredes 2005 ), although some authors have 
proposed that many if not all REXBs could be MQs (e.g. Fender 
2004). Thus, the real number of MQs in our galaxy could be 
notably higher. 

MQs are considered scaled-down versions of distant quasars 
since they present many of the characteristics of these extra- 
galactic sources. They serve as suitable scenarios for understand- 
ing a number of processes, such as mass accretion onto the com- 
pact object or jet formation and evolution, in timescales inacces- 
sible in the case of quasars. In analogy with radio-loud quasars, 
which show jets impacting on the intracluster medium, one may 
expect to find radiative signatures due to strong shocks in the 
termination regions of MQ jets as well. Hot spots and double- 
lobe morphologies are common features of the powerful extra- 



galactic FRII sources (Faranoff & Riley 119741 ). in which a va- 
riety of large-scale non-thermal structures are revealed at radio 
wavelengths. In the case of MQs, convincing evidence of jet in- 
teraction with the interstellar medium (ISM) are not numerous, 
partially due to the relatively small number of known sources. 
Nevertheless, hints or evidence of such interactions have been 
observed, at different wavelengths, in SS 433 (Zealey et al. 
[19801 IE 17 40.7- 2942 (Mirabel et al. [19921) , XTE J1550-564 
(Corbel et al.| 2QTj2|, Cygnus X-3 (Hein dl et al. 120031 , Cygnus X- 
1 (M artlet a l. [19961 Gallo et al. 120051 ), GRS 19 15+105 (Kaiser 
et al. 120041 ). HI 743 - 322 (Corbel et al. 120051 ). LS I + 61 303 
(Paredes et al. I2007afa ). and Circinus X-l (Tudose et al.[2006 ), 
although only some theoretical work has been done regard- 
ing the hydrodynamics of the interaction, or the resulting non- 
thermal emission (e.g. Aharonian & Atovan fl9981 Velazquez & 
Raga 2000, Heinz & Sun yaev 120021 Bosch-Ramon et al. 120051 
Perucho & Bosch-Ramon |2008l ). 

MQ ejections can transport large amounts of kinetic energy and 
momentum very far from the binary system. The age of the 
source times the jet kinetic power, tmq X Qjet, can be as high 
as 10 12 s (3 x 10 4 yr) x 10 37 - 10 3 9 erg s" 1 ~ 10 49 - 10 51 erg 
(see e.g., Cygnus X-l, Gallo et al. 120031 SS 433, Zealey et al. 
1980). Even for low radiative efficiencies converting this energy 
into non-thermal emission, MQ jet termination regions could 
produce significant fluxes if they were at distances of few kpc. 



2 



P. Bordas, V. Bosch-Ramon, J. M. Paredes and M. Perucho: Non-thermal emission from microquasar/ISM interaction 



It could well be that MQs were associated, via interaction with 
their surrounding medium, to some of the extended non-thermal 
radio sources detected in the Galaxy of unknown origin (e.g., 
Paredes et al. 12007V ). 

In this work, we investigate whether a typical MQ fulfills the 
conditions to be detectable when interacting with the surround- 
ing external gas. We adopt some reasonable assumptions for the 
MQ power and age, and the density of the ISM, as well as for the 
non-thermal energy and magnetic field equipartition fraction in 
the interaction regions. We study the case of a high-mass system 
in order to see the role of a massive and hot companion. In the 
case of a low-mass system, the radiation field of the companion 
star and therefore the Inverse Compton (IC) contribution would 
be much lower. 

This paper is organized as follows. In Section 2, we present an 
analytical interaction model to characterize the three shocked 
zones; for each region, the main properties of the non-thermal 
emitters are studied, and the physical conditions and geometry 
are established. In Sect. 3, we calculate the total emission pro- 
duced through synchrotron, relativistic Bremsstrahlung, and IC 
processes for the different set of parameters used in the model. 
We carried out hydrodynamical simulations to study in more de- 
tail the jet evolution when interacting with the ambient medium. 
These simulations allow us also to test and enforce the validity of 
the assumptions of the analytical model. This is done in Sect. 4. 
The details of the hydrodynamical simulations can be found in 
Appendix A. Finally, in Sect. 5 we discuss the obtained results, 
giving detectability predictions at different energy bands. 

2. Model description 

We use a model adapted from those applied to extragalactic 
sources in order to investigate the interactions of MQ jets with 
their environment. Falle ( 1991 ) and Kaiser & Alexander ( 1 19971 1 
have developed an analytical model where self-similarity is used 
to characterize the evolution of the jet/medium interaction struc- 
tures in radio galaxies. We have closely followed their work, 
with minor modifications required for the application of the 
model to the MQ context. 

We consider twin conical jets emerging from the central source 
in opposite directions. The ejecta begin to decelerate when the 
mass of the swept up external gas becomes similar to the mass 
carried by the jet. Two shocks are formed at the jet tip: a forward 
shock (the bow shock) propagating into the ISM and a reverse 
shock directed inwards into the jet material. Matter crossing the 
reverse shock inflates the cocoon, which protects the jet from 
disruption due to turbulent gas entrainment. Moreover, a recon- 
finement shock is also formed in the jet at the point where its 
pressure equals that of the cocoon. 

In our scenario, all the shocks are assumed to be strong. The 
density and pressure of the shocked material are calculated us- 
ing the Rankine-Hugoniot conditions for a strong shock at each 
front. We note that since the conditions in the shocked regions 
evolve with time, all related physical variables will also change 
with time. Therefore, the non-thermal particle population is cal- 
culated considering the time dependence of the radiative and adi- 
abatic losses (see Sect. 2.3). Given that all the considered shocks 
are weakly relativistic, we adopt an adiabatic index y = 5/3. 
We assume the presence of a randomly oriented magnetic field 
B in the downstream regions, derived taking the magnetic en- 
ergy density to be ~ 10 % (77 = 0.1) of the thermal energy den- 
sity. In each shock region, the fraction of kinetic luminosity as- 
sumed to be transferred to non-thermal particles is taken to be 
1 % (x = 0.01). 
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Fig. 1. Scheme of the model (not to scale) representing the three 
different zones. The central binary system is located at the cen- 
ter. Two jets emerge from it and extent outwards until they are 
effectively decelerated at a distance ~ 10 20 cm. Jet material that 
crosses the reverse shock inflates the cocoon, which expands ex- 
erting work on the shocked ISM. A contact discontinuity sepa- 
rates the cocoon and the shocked ISM, the latter being further 
limited by the bow shock. 



2. 1 . The properties of the three emitting zones 

We focus our studies in three separated emitting zones. The first 
one, the shell region, corresponds to the ISM material swept up 
by the bow shock. The second one, the cocoon region, corre- 
sponds to the material of the jet crossing the reverse shock. The 
third zone accounts for the shocked jet material after crossing the 
conical reconfinement shock. A sketch of the model is shown in 
Fig.[TJ 



2.1 .1 . Bow shock 

We are interested in the case when the bow shock can accel- 
erate particles, and we concentrate on its adiabatic phase (see 
Sedov [T959] ), which implies that the inertia of the medium must 
be large enough to take a significant fraction of the energy and 
momentum from the jet. Although in the extragalactic context 
it is thought that the non-thermal radio emission comes solely 
from the shocked jet material in the cocoon, it is worth to explore 
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the possibility to have also a contribution from the shocked ISM 
material. After the interaction structures have reached a charac- 
teristic length-scale / = [Q 2 Jp 2 SM c 6 (T iet - 1) 3 ] 1/4 , where Q n 
is the jet power, pism is the ISM mass density, Lj e t is the jet bulk 
Lorentz factor and c is the speed of light, the equations for the 
length and velocity of the bow shock at a given source age ?mq 
are (see Falle ( fTWO , Kaiser & Alexander dl997l ): 



Table 1. List of the parameters which remain with a constant 
value in the analytic model 



k = c 
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respectively. c\ ~ 1 is a dimensionless constant that depends 
on the adiabatic index of the jet material and the geometry of 
the bow shock. The pressure in the shell can be calculated as 

Pb = (3/4)pi SM v 2 h . 

To determine the radius of the bow shock, rb, we assume a 
self-similar relationship between its length and width given by 
R = Zb/n>, which can change from source to source in the range 
R e [2 - 3 ] (Leahy et al. [19891 Kaiser et al. [20041 Perucho & 
Marti 2007). Here we use a fiducial value of R = 3 (similar to the 
value from our simulations; see Appendix |AV We assume here 
a plane shock and do not enter in details concerning the precise 
geometry of the bow shock. Nevertheless, we need the accelera- 
tor size in order to compute the maximum energies that particles 
can reach, as well as an estimation of the emitter size to provide 
the surface brightness of the source. In the case of the bow shock, 
both the accelerator and the emitter size are taken to be ~ r/,. The 
mass density and the pressure, which has been assumed equal to 
that of the cocoon, have been taken homogeneous in the whole 
shell. Since the bow shock velocity evolves with time, the mass 
density and the pressure change as well. Therefore, the magnetic 
energy density varies also with time. Finally, the radiation field 
energy density depends on u+ = L*/4 n l^c (where is the 
companion luminosity). 

2.1 .2. Cocoon 

The cocoon is filled with jet material that crosses the reverse 
shock. A contact discontinuity separates the cocoon and the 
shell region materials. We do not consider the mixing of the gas 
of both zones that could occur due to Kevin Helmholtz and/or 
Rayleigh Taylor instabilities. The cocoon pressure is taken to be 
the same as that in the other side of the contact discontinuity at 
any time, P c = Pb. Given the strong compression of the ISM gas 
in shell zone, the length and the width of the cocoon are taken 
to be also ~ 4 and ~ r\,, respectively. For the accelerator size, 
the reverse shock, we adopt the constant jet width after the re- 
confinement point (see below). The pressure and mass density 
of the material can be taken as homogeneous for the whole co- 
coon to first order approximation. The homogeneity assumption 
comes from the fact that the cocoon material is expected to be 
subsonic. Both pressure and mass density will depend on the 
time evolution of the reverse shock velocity, like the magnetic 
and radiation fields, in the same way as in the case of the shell. 



2.1.3. Reconfinement region 

The conical jet initially emerging from the central engine has 
an opening angle ¥ =s r/z ~ 0.1 rad, where z is the distance 



Parameter 


Symbol 


Value 


Jet Lorentz factor 


Ljet 


1.25 


Jet half opening angle 




0.1 rad 


Luminosity companion star 




10 39 erg s - 1 


Self-Similar parameter 


R 


3 


Magnetic equipartition fraction 


'1 


0.1 


Non-thermal luminosity fraction 


X 


0.01 



to the injection point and r is the jet radius. Given the pressure 
exerted by the surrounding cocoon, the jet radius becomes ap- 
proximately constant at z CO nf , the distance where the lateral jet 
pressure becomes similar to the cocoon one. Following Kaiser 
& Alexander (1997), we obtain the following radius for the re- 
confined jet: 



2Qjet v jet 



\(y+l)(r jet -l)7r C 2 P c ' 



(3) 



where y — 5/3 is the adiabatic index of the cocoon material and 
vjet is the jet velocity. For a recent semi-analytical treatment of 
reconfinement shocks, see (2008). After injection downstream, 
the shocked material still moves almost at the jet velocity. The 
normal component of the velocity of the reconfinement shock, 
i.e., the one that suffers a discontinuity, is much smaller than vj et . 
The jet keeps most of its thrust and remains supersonic until it 
is braked at the reverse shock. From the reconfinement point the 
jet roughly keeps a constant radius. We neglect further recolli- 
mation shocks that may occur (see AppendixlAl. 
The size of the accelerator is taken as the width of the jet at 
the recollimation shock, r con f . Once accelerated, relativistic par- 
ticles are approximately advected at the jet velocity, and prop- 
agate up to the reverse shock. The length of the emitter, as- 
sumed to be one-dimensional, is therefore taken to be ~ l\, since 
Zconf ~ '"conf/^ <k k- The magnetic field in the reconfinement 
region is derived as in the other interaction zones. The radiation 
density, however, gradually decays from z con f to lb as the inverse 
square of the distance to the companion star. The density and 
the pressure are considered to remain roughly constant after the 
reconfinement point since the jet radius is constant beyond z con f • 



2.2. The non-thermal particle populations 

The theory of diffusive shock acceleration in the llinear regime 
for non-relativistic velocities predicts a power-law index of p ~ 
2 for the injected particle spectrum in strong shocks (e.g. Drury 
119831 ). We take this value for the non-thermal particle spectra 
injected at the reconfinement, reverse and bow shock fronts, i.e. 
N(E) = KE 2 . As noted above, the normalization constant K is 
taken such that ~ 1 % of the kinetic power flowing in the jets is 
converted into non-thermal energy in the postshock region right 
after each shock. We remark that the radiation luminosities scale 
linearly with this non-thermal fraction. 

The maximum energies of the relativistic particles, E m . dx , are cal- 
culated at any source age equating the energy gain to the differ- 
ent cooling processes. In case radiative cooling is not effective, 
the maximum energy is constrained by the Hillas limit (Hillas 
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1984), i.e. the particle gyroradius, r g , equals the accelerator size. 
We have adopted Bohm diffusion (D - r g c/3) and the magnetic 
field in each shocked region, to calculate the rate at which parti- 
cles gain energy (e.g. Protheroe 1 19991 1: 

£accei(0 « (3/20) [v s (t)/c] 2 e B(t) c . (4) 

Focusing hereafter on electrons, the particle energy distribution 
is computed taking into account radiative (synchrotron, rela- 
tivistic Bremsstrahlung and IC emission; see, e.g., Blumenthal 
& Gould 1970) and adiabatic losses along the lifetime of the 
source. As pointed out in Blundell et al. (2000) in their study 
on extragalactic double radio sources, the spectral aging of the 
non-thermal populations depends on the evolution of the physi- 
cal conditions in each region, in particular the expansion veloci- 
ties and magnetic fields for each interaction zone. At each time, 
the conditions for the radiative and adiabatic losses change. The 
non-thermal particle distribution at a given source age ?mq has to 
be therefore computed summing up the differently evolved con- 
tributions of the particles injected previously at the shock fronts 
all along the source lifetime until reaching fMQ- 
For the synchrotron losses, we use the magnetic field considered 
above for each interaction region. Relativistic Bremsstrahlung 
is calculated accounting for the densities in the downstream 
regions. To compute IC losses, we consider the dominant ra- 
diation field from the companion, an OB star with luminosity 
L* = 10 39 erg s _1 . Adiabatic losses, E (t) - [v(f) / r (t)] E, are 
computed from the size of the emitters r(f) and the expansion 
velocity v(f). At the downstream region, after the reconfinement 
shock, there is no expansion since the jet radius keeps roughly 
constant. In the case of the shell region, the expansion velocity is 
Vb, and for the cocoon, it is ~ 3/4 Vb (Landau & Lifshitz 1987). 
As noted above, the electron population properties, the adiabatic 
coefficient, and the magnetic and radiation fields are taken ho- 
mogeneous as a first order approximation in the bow shock and 
cocoon regions. A study of the detailed spatial structure of the 
magnetic field, pressure and mass density at each interaction re- 
gion is beyond the scope of this work. For the electrons injected 
at the reconfinement shock, we compute their evolution assum- 
ing a one-dimensional emitter, the jet, with the stellar radiation 
density decaying as oc l/z 2 . 

The list of the parameters which remain constant in our analytic 
model are pesented in Table 1 . 

Given £ mal and the evolved electron distribution in each emitting 
region, and accounting for the mentioned radiation mechanisms, 
we have computed the SEDs for each zone. This is presented in 
the next section. 

3. Model Results 

We have studied how the computed SEDs are affected by vary- 
ing the source age, ?mq, the J et kinetic power, Q jet , and the ISM 
density, «ism- In Figs. [2] to [5] we show the SEDs for the shell 
(top), cocoon (middle) and jet reconfinement (bottom) regions. 
The contribution of only one jet interacting with the ISM is ac- 
counted. The whole set of parameter values, required for the 
different cases to obtain the particle evolution and subsequent 
emission, are presented in Table 2. 

Synchrotron emission from the three interaction zones is the 
channel through which the highest radiation output is ob- 
tained, with bolometric luminosities that can be as high as 
~ 10 33 erg s _1 for powerful sources (see Figs. 4 and 5). At 
high and very high energies, IC emission is the dominant pro- 
cess in the cocoon and reconfinement regions, reaching lumi- 
nosities up to ~ 10 30 erg s , while for the shell zone relativistic 



Bremsstrahlung dominates at this energy range, with luminosi- 
ties up to ~ 10 32 erg s -1 . Notable differences are found in the 
reported SEDs when varying the source age ?mq from 10 4 to 
10 5 yr. For older sources, the interaction zones are located at 
larger distances from the companion star, its radiation energy 
density m+ decreases and the IC contribution gets slightly lower, 
although particle aging leads to higher emission at lower ener- 
gies in the shell and cocoon regions. Relativistic Bremsstrahlung 
in the shell zone is strongly affected by the source age, due to 
a particle accumulation effect, giving luminosities a factor of 
~ 10 larger for old sources. Regarding the ambient medium den- 
sity, higher values of «ism make the jet to be braked at shorter 
distances from the central engine. The interaction regions are 
then under a higher radiation energy density from the compan- 
ion star, and the IC emission is accordingly enhanced. The rela- 
tivistic Bremsstrahlung emission in the shell zone is also higher 
for denser mediums, since the luminosity is proportional to the 
target ion field density, n t = 4«ism- Finally, we note that in our 
model all the non-thermal luminosities are proportional to 2j et . 

3. 1 . Emission from the shell region 

Under the acceleration conditions in this region, the maximum 
energies £ ma , of electrons for the different cases range 2-10 TeV, 
being always limited by synchrotron losses. Synchrotron emis- 
sion is also the dominant radiative channel, peaking at higher 
frequencies for higher E max . The electron energy distribution has 
an energy break, E\» at which the radiative cooling timescale and 
the age of the source are the same. E^ changes its value depend- 
ing on the physical conditions: older sources render smaller E^, 
since particles have cooled via synchrotron for a longer time; 
higher densities yield lower E^, due to smaller 1^ and therefore 
larger B and synchrotron cooling. Larger jet powers enhances 
B as well, reducing again E\, (this discussion about E^ is valid 
also for the other interaction regions). The energy break is man- 
ifested in the synchrotron spectrum at frequencies in the range 
~ 10 14 — 10 16 Hz. Given the long synchrotron cooling timescales 
T sync ~ lO 14 ^:^ 1 ^ s (B_ 4 = fi/10~ 4 G) compared with the 
age of the source ~ 10 12 s, only the highest energy band of the 
particle spectrum has reached the steady state. The synchrotron 
luminosities reach values up to ~ 10 33 erg s _! (see Figs. |4] and 
0. The Relativistic Bremsstrahlung has a strong dependence on 
the density, follows the energy distribution of electrons, and can 
reach luminosities ~ 10 32 erg s _1 in the 1 MeV-10 GeV range 
(see Fig. O. The IC luminosities are well below the synchrotron 
ones in all the studied cases . 

3.2. Emission from the cocoon 

Electrons here can reach energies of about E mlix ~ 100 TeV, be- 
ing limited by escape losses. Synchrotron emission, again the 
dominant radiative channel, extends up to hard X-rays in the co- 
coon, with luminosities reaching up to ~ 10 33 erg s _1 (see Fig. I5J 
depending on the different parameter values. As in the shell re- 
gion, the IC luminosities are well below the synchrotron ones. 
Otherwise, due to the low densities in the cocoon, relativistic 
Bremsstrahlung is negligible in this zone. 

3.3. Emission from the jet reconfinement region 

The size of the accelerator (Hillas 1984) put a limit of ~ 10 TeV 
on the maximum energies that electrons can acquire at the re- 
confinement shock. An increase in both the synchrotron and the 
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Table 2. Parameter values adopted to compute the SEDs for the three emitting zones. The target density, n t , is only shown when 
non-negligible. We show the parameter values at f M Q, but we remark that they smoothly vary with time. The computation of the 
non-thermal particle distribution have taken this into account. See Section 2 for details. 



Parameter 


Fig. 2 


Fig 


;.3 


Fig 


;.4 


Fig. 5 


Jet kinetic power g jct (erg s "' ) 


10 36 


10 37 


ISM density /2ism (cm - "*) 


0.1 


1 


0.1 


1 


Source age f M o (yr) 


10 4 


10 5 


10 4 


10 5 


10 4 


10 5 


10 4 


10 5 


SHELL 




Magnetic field B (G) 


2.9 xlO" 5 


1.1 xlO" 5 


5.8 xlO" 5 


2.3 xlO" 5 


4.6 xlO" 5 


1.8 xlO" 5 


9.2 xlO" 5 


3.6 xlO" 5 


Shock velocity v b (cm s" 1 ) 


4.3 x 10 7 


1.7 x 10 7 


2.7 x 10 7 


1.0 x 10 7 


6.9 x 10 7 


2.7 x 10 7 


4.3 x 10 7 


1.7 x 10 7 


Emitter size r (cm) 


2.3 xlO 19 


9.1 xlO 19 


1.4 xlO 19 


5.7 xlO 19 


3.62 xlO 19 


1.44 xlO 20 


2.3 xlO 19 


9.1 xlO 19 


Rad pnpr^v dens u . (pyp cm"^ 


5.0 xlO" 12 


3.2 xlO" 11 


1.2 xlO" 13 


8.0 xlO" 13 


2.0 xlO" 12 


1.3 xlO" 12 


5.0 xlO" 13 


3.2 xlO" 13 


Maximum enerpv F fTe\H 


8.1 


5.1 


3.6 


2.3 


10.2 


6.4 


4.5 


2.8 


Target density n t (cm~^) 


0.4 


0.4 


4.0 


4.0 


0.4 


0.4 


4.0 


4.0 


COCOON 




Magnetic field B (G) 


3.8 xlO" 5 


1.5 xlO" 5 


7.5 xlO" 5 


3.0xl0- 5 


5.9 xlO" 5 


2.4 xlO" 5 


1.2 xlO" 4 


4.7 xlO" 5 


Shock velocity v s (cm s" 1 ) 


1.8 xlO 10 


1.8 xlO 10 


1.8 xlO 10 


1.8 xlO 10 


1.8 xlO 10 


1.8 xlO 10 


1.8 xlO 10 


1.8 xlO 10 


Emitter size r (cm) 


6.5 xlO 18 


2.5 xlO 18 


4.0 xlO 19 


1.6 xlO 19 


1.0 xlO 19 


4.0 xlO 18 


6.5 xlO 19 


2.5 xlO 19 


Rad. energy dens. (erg cm" 3 ) 


7.0 xlO" 12 


1.7 xlO" 11 


4.4 xlO" 13 


1.1 xlO" 12 


2.8 xlO" 12 


1.7 xlO" 12 


7.0 xlO" 13 


4.4 xlO" 13 


Maximum energy £ max (TeV) 


275.5 


275.5 


275.5 


275.5 


871.2 


871.2 


871.23 


871.2 


RECONFINEMENT 




Magnetic field B (G) 


2.6 xlO" 5 


1.0 xlO" 5 


5.2 xlO" 5 


2.1 xlO" 5 


4.2 xlO" 5 


1.6 xlO" 5 


8.3 xlO" 5 


3.3 xlO" 5 


Shock velocity v conf (cm s" 1 ) 


1.8 xlO 9 


1.8 xlO 9 


1.8 xlO 9 


1.8 xlO 9 


1.8 xlO 9 


1.8 xlO 9 


1.8 xlO 9 


1.8 xlO 9 


Emitter size r (cm) 


1.9 xlO 19 


7.6 xlO 19 


1.2 xlO 19 


4.8 xlO 19 


3.0 xlO 19 


1.2 xlO 20 


1.9 xlO 19 


7.6 xlO 19 


Rad. energy dens, (erg cm" 3 ) 


2.5 xlO" 8 


4.0 xlO" 9 


1.0 xlO" 7 


1.6 xlO" 8 


6.4 xlO" 9 


1.0 xlO" 9 


2.5 xlO" 8 


4.0 xlO" 9 


Maximum energy £ max (TeV) 


3.2 


3.2 


3.2 


3.2 


10.1 


10.1 


10.1 


10.1 



IC emissivities is found for denser mediums, in which case P c 
is larger yielding smaller z con f, which determines the magnetic 
and radiation field energy densities. The synchrotron luminosity 
is similar to the IC one for weak jets, whereas it is well above it 
for powerful jets, where luminosities ~ 10 32 erg s" 1 are obtained 
(see Figs. 4 and 5). This is due to the stronger dependence of 
the synchrotron emission with the jet power. The stellar radia- 
tion field enhances the IC contribution when the set of parame- 
ters make z con f shorter. As in the cocoon region, the relativistic 
Bremsstrahlung contribution is negligible. 



4. Hydrodynamical simulations vs analytical model 

We performed hydrodynamical simulations to further study the 
interaction of MQ jets with their surroundings (see Appendix 
for details). The values of Qj et , ?mq and «tsm adopted in the hy- 
drodynamical simulations are similar to those assumed in the 



analytical model, so it is worthwhile checking the main phys- 
ical quantities used to compute the non-thermal emission. The 
pressure evolution with simulation time have been computed for 
the shell and cocoon regions. We show also a fit of their evolu- 
tion performed to check their values up to ?mq — 10 5 yr. From 
an initial value of ~ 2 x 10" 9 erg cm" 3 , the cocoon pressure 
smoothly reaches ~ 8 x 10" 11 erg cm" 3 at t = 7 x 10 11 s. 
The shell behaves similarly, as expected, with a final pressure of 
~ 8 x 10" 11 erg cm" 3 . These values are in reasonable agree- 
ment with those found in the analytical model, ranging between 
(2 - 7) x 10" 10 erg cm -3 . Concerning the mass densities, pism 
stabilizes at ~ 2 x 10" 25 g cm" 3 and p co ~ 4 x 10" 29 g cm" 3 
for the shell and the cocoon regions, respectively. These values 
are again very similar to those found in the analytical treatment. 
On the other hand, the geometry of the emitting structures is 
related to the self-similar ratio R. In the analytical model, we 
use R = 3. We find this value to be in accordance with the re- 
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Fig. 2. Obtained SEDs for the shell (top), cocoon (middle) and 
jet reconfinement (bottom) regions taking the values for the jet 
power (2jet = 10 36 erg s _! and an external gas particle density 
"ism = 0.1 cm 4 . Two different values for the source age are 
represented, ?mq = 10 4 yr (solid line) and 10 5 yr (dashed line). 
See Table 2 for details on the parameter values for each region. 



suits of the numerical simulations, which yield a value between 
2.5 and 2.7. Finally, we note that in our model only a strong 
shock at the reconfinement point is assumed, and no additional 
shocks are considered along the jet. Otherwise, the hydrodynam- 
ical simulations show the existence of several conical shocks that 
develop in the jet when its pressure falls to that of the surround- 
ing cocoon (see Fig. |A. lb . Therefore, the non-thermal emission 
presented in Figs.|2]to|5]for the reconfinement region, with only 
one strong shock and an acceleration efficiency of x — 0.01, 
should be taken as a rough approximation of the real situation. 



Fig. 3. Same as in Fig. 2 but taking «ism = 1 cm 3 . See the 
physical parameter values in Table 2 



5. Discussion 

We have explored whether non-thermal emission can be ex- 
pected from the bow shock, the cocoon and the reconfinement 
regions in the interaction of MQ jets with their environment. 
In the extragalactic framework, non-thermal radiation is sup- 
posed to come only from the cocoon as extended radio emis- 
sion, as well as locally in the hot spots at the jet tips. Although 
the shell of shocked ambient material plays an important role in 
the analytical models describing the growing of FR II galaxies 
(Blandford et al. 1 19741 Scheuer fl974| l, no radio emission is usu- 
ally assumed to be produced there (see, however, Rudnick et al. 
1980). In our model, the bow shock velocity is still large enough 
to accelerate particles, so it is worth to compute the expected 
non-thermal emission also from this region. Our results show 
a contribution from this zone that is comparable to that of the 
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Fig. 4. Same as in Fig. 2 but taking Qj et = 10 erg s . See the Fig. 5. Same as in Fig. 2 but taking Qjet = 10 erg s and 



physical parameter values in Table 2 



«ism = 1 cm . See physical parameter values in Table 2 



cocoon and higher than that coming from the reconfinement re- 
gion. The geometry of the interaction structures, however, makes 
it difficult to disentangle the emission from the cocoon and the 
bow shock region since they could appear co-spatial in the plane 
of the sky. A clear relativistic Bremsstrahlung would favour a 
shell origin of the emission. 

The highest radiation output within the studied set of parame- 
ters corresponds to the case: Qj et = 10 37 erg s _1 , ?mq = 10 5 yr 
and «ism = 1 cirT 3 . In case the emitting source were located 
at 3 kpc, this would imply a flux density of ~ 150 mJy at 
5 GHz. The emitting size would be of a few arcminutes, since 
the electron cooling timescale is longer than the source lifetime 
and they can fill the whole cocoon/shell structures. Considering 
this angular extension and taking a radio telescope beam size 
of 10", radio emission at a level of ~ 1 mJy beam -1 could be 
expected. At the X-ray band, we find a bolometric flux in the 



range 1-10 keV of i*ViokeV ~ 2 x 10~ 13 erg s~' crrT 2 . The elec- 
trons emitting at X-rays by synchrotron have very short time- 
scales, and the emitter size cannot be significantly larger than 
the accelerator itself. Although the X-rays produced in the shell 
through relativistic Bremsstrahlung are expected to be quite di- 
luted, the X-rays from the cocoon would come from a relatively 
small region close to the reverse shock, and could be detectable 
by XMM-Newton and Chandra at scales of few arcseconds. In 
the gamma-ray domain, the flux between 100 MeV and 100 
GeV is F 100 MeV<E<ioo Gev ~ 10~ 14 erg s~' crrT 2 , while the inte- 
grated flux above 100 GeV is FE>iooGeV ~ 10~ 15 erg s~' cm -2 . 
These values are too low to be detectable by current Cherenkov 
telescopes. Taking into account the rough linearity between 
Gjet, «ism, ty[Q, x and dT 1 with the gamma-ray fluxes obtained, 
sources with higher values of these quantities than the ones used 
here may render the MQ jet termination regions detectable by 



8 



P. Bordas, V. Bosch-Ramon, J. M. Paredes and M. Perucho: Non-thermal emission from microquasar/ISM interaction 



present Cherenkov telescopes. For the weakest jets, i.e, low- 
est ISM densities and youngest sources adopted in our model 
(Gjet = 10 36 erg s _1 , f M Q = 10 4 yr and n ISM = 0.1 cm- 3 ), the 
fluxes are strongly suppressed. In the radio band, the specific 
flux is F5 ghz ~ 0.1 mJy beam -1 ; the integrated flux at X-rays 
in this case is fViokeV ~ 10~ 14 erg s _1 crrT 2 , and at gamma- 
rays F iooMeV<E<iooceV ~ 2 x 10~ 17 erg s" 1 cm" 2 and F E >iooGeV ~ 
2 x 10~ 16 erg s" 1 cm" 2 . 

We remark that the fluxes showed above strongly depend on the 
non-thermal luminosity fraction x- In the present model we use 
a quite conservative value x - 0.01. In the case of a source 
able to accelerate particles at a higher efficiency at the interac- 
tion shock fronts, then the expected non-thermal fluxes would be 
enhanced by a factor (x/0.01). 

The non-thermal emission from the interaction regions in MQs 
presents characteristic features that can be distinguished from 
emission coming from the central binary system. First of all, the 
interaction structures are localized at distances up to ~ 10 20 cm. 
The lifetime of electrons radiating at high and very high ener- 
gies is relatively short thus the emission region should be lo- 
calized near the accelerator (but far away from the central sys- 
tem). Regarding radio emission, synchrotron cooling times are 
expected to be larger than the source lifetime, r MQ = 10 4 - 10 5 yr. 
The cocoon and bow shock emitting regions would be expected 
to form a kind of radio nebula around the central system and the 
flux densities at the level showed above would then come from 
a quite extended source. 

The detection of non-thermal emission from the interaction 
zones would be a proof that efficient acceleration of particles up 
to relativistic energies is taking place far away from the central 
binary system. The averaged kinetic power carried away in the 
jets could be better constrained, eventually showing that it can 
be much higher than that inferred directly from observations of 
the inner jet emission alone (Gallo et al. 2005) Heinz [2002 ) 
Despite we focus on the non-thermal emission from the MQ jet 
termination regions, thermal Bremsstrahlung emission should 
be expected from the shell. Although the shocks considered 
here are still adiabatic, a non-negligible fraction of the jet ki- 
netic luminosity of up to a few % may be radiated via ther- 
mal Bremsstrahlung. For bow-shock velocities of few times 
10 7 cm s _1 , the thermal emission would peak at UV-soft X-rays, 
energy band that is strongly affected by absorption in the ISM. 
Observations of the thermal radiation in radio and optical from 
the interaction structures have been used to extract information 
of the shell physical conditions (e.g. Cygnus X-l, Gallo et al. 
120051 Russell et al. l2007l . 

The role of accelerated protons in the shocks deserves a few 
words, since it may be relevant in some specific situations. Given 
the conditions in the strong shocks we are considering, relativis- 
tic protons may reach energies of about 100 TeV; for shell den- 
sities ~ 1 cirT 3 , the accelerated protons have lifetimes of about 
10 15 s. To reach the gamma-ray fluxes detectable for the present 
generation of satellite borne and ground based gamma-ray in- 
struments above ~ 100 GeV, ~ 10" 13 erg cirT 2 s _1 , the energy 
in non-thermal protons stored in the shell should be as high as 
~ 10 48 erg at few kpc distances. For a source age of ?mq = 10 5 yr, 
the required injected power in relativistic protons should be 
about ~ 3 x 10 35 erg s _1 , thus implying that moderate levels 
of hadronic emission from MQ jet termination regions may be 
eventually detected from powerful jets, i.e. Q ]et Z 10 37 erg s . 
The reason why some MQs show non-thermal emission from the 
jet/ISM interaction regions, whereas in other cases such emis- 
sion remains undetected, is still unclear. In the context of our 
model, we can study the effects of varying the set of parame- 



ters defining the source and their environment, and predict some 
cases in which the interaction structures may or may not be de- 
tectable. First of all, the energy input injected to the medium 
should be high enough, and there exist strong differences in the 
jet kinetic power from source to source. In addition, it could 
be also the case that the density of the surrounding medium is 
so low that the shell and the cocoon get very large and their 
radiation too diluted to be detectable (see, for instance, Heinz 
2002). Moreover, MQ jets could get disrupted at some source 
age, as it is found in FR I galaxies. If this happened within times 
much shorter than the MQ lifetime, the probability to detect a 
cocoon/shell structure in the MQ surroundings would be smaller 
(although the detection of some other kind of structures is not 
discarded). On the other hand, some sources may be too far, or 
the non-thermal fraction too small, to detect significant emission 
from the interaction regions. 

The evolution of the pressure, mass density, the velocities and 
the Mach number predicted by the analytical model are in good 
agreement with those found in the hydrodynamical simulations 
for the shell and the cocoon regions. Otherwise, these simula- 
tions show that several conical shocks may be present within the 
jet as a consequence of pressure adjustments with the surround- 
ing cocoon, instead of the one strong shock adopted in the ana- 
lytical treatment. Finally, the length and width of the structures 
in the model and those found through the numerical simulations 
are also similar, with a constant ratio R ~ 3 in both cases, imply- 
ing that the physical assumptions used in the analytical treatment 
are valid to first order. 

The results of this work show that the surroundings of some MQs 
could be extended non-thermal emitters from radio to gamma- 
rays. In addition, from a comparison with observations, the mag- 
netic field and the particle acceleration efficiency in the jet/ISM 
interaction regions can be constrained, giving an insight on the 
physics of these interaction structures. It is interesting to note 
that, although the adopted model is rather simple, it already ac- 
counts for cases when the sources should remain undetectable 
and cases in which radiation could be detected. 
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Appendix A: hydrodynamical simulations 

The simulation has been performed in order to check the phys- 
ical values adopted in the analytical model. We used a two- 
dimensional finite-difference code based on a high-resolution 
shock-capturing scheme, which solves the equations of relativis- 
tic hydrodynamics written in conservation form. This code is an 
upgrade of the code described in Marti et al. ( 1997) (see Perucho 
et al. l20051 >. The simulation has been carried out in two dual-core 
processors in the Max-Planck-Institut filr Radioastronomie. 

The numerical grid of the simulation is formed by 320 cells 
in the radial direction and 2400 cells in the axial direction in a 
uniform region, with physical dimensions of 40 x 600 rj. An ex- 
panded grid with 320 cells in the transversal direction, brings 
the boundary from 40 rj to 500 rj, whereas an extended grid in 
the axial direction, consisting of 440 extra cells, spans the grid 
axially from 600 rj to 900 r y The enlargement of the grid is nec- 
essary to ensure that the boundary conditions are sufficiently far 
from the region of interest, and to avoid numerical reflection of 
waves in the boundaries affecting our results. The conditions at 
the boundaries are reflection on the jet axis and in the side where 
the jet is injected, simulating the presence of the counter-jet co- 
coon, with the exception of the injection point, where inflow 
conditions are used. Finally, outflow conditions in the outer axial 
and radial boundaries are used. 

The numerical resolution in the uniform grid is thus of 8 
cells/rj in the radial direction and 4 cells/rj in the axial direc- 
tion. The low resolution used is justified by the fact that we are 
interested in the macroscopic features of the jet and backflow, 
but not in mixing and turbulence, allowing much less time con- 
suming simulations. All the physical variables are scaled to the 
units of the code, which are the jet radius r, the rest-mass density 
of the ambient medium pism, and the speed of light c. 

The jet is injected in the grid at a distance of 10 18 cm from the 
compact object, and its initial radius is taken to be ro = 10 17 cm. 
The time unit of the code is thus equivalent to « 3 x 10 6 s, which 
is derived using the radius of the jet at injection and the speed 
of light (ro/c). Both the jet and the ambient medium are con- 
sidered to be formed by a non-relativistic gas with adiabatic ex- 
ponent T = 5/3. The number density in the ambient medium is 
"ism = 0.3 ctrT 3 and the temperature is T — 100 K. The role of 
the temperature is not relevant in our context. The velocity of the 
jet at injection is 0.6 c, its number density n s = 1.4 x 10" 5 cm" 3 , 
and temperature T ~ 10 11 K (which corresponds to a sound 
speed ~ 0.1 Vj e t)- These parameter values result in a jet power 
gje t = 3 x 10 36 erg s _1 . Figs. IA.4I and IA.3I show the velocity, 
Mach number, mass density and pressure maps obtained with 
the numerical simulations. The upper (lower) panels in Fig. IA.5I 
show the evolution of the pressure and the mass density with 
simulation time for the shell (cocoon) region. The evolution of 
the self-similar parameter R reaches a value between 2.5 and 2.7 
at the end of the simulations, as can be seen in Fig. IA.2l 

At the time the simulation is stopped, after evolving as 2.3 x 
10 4 yr, the bow-shock is moving at a speed » 3 x 10 7 cms _1 , 
and has reached a distance ~ 3.3 x 10 19 cm. Initially, the jet ex- 
pands, accelerating and cooling, due to its initial overpressure. 
When the flow becomes underpressured with respect to the am- 
bient medium, a first reconfinement shock is generated close to 
the injection 2 x 10 1 8 cm. The fluid becomes then again over- 
pressured when passing through the shock and this process is 



repeated several times around pressure equilibrium with the ex- 
ternal medium. The subsequent shocks produced by these oscil- 
lations around equilibrium are stronger than the first, suggesting 
some coupling to a pinching Kelvin-Helmholtz instability. The 
jet is supersonic at injection, with Mach number Mj = 6.5 on 
axis, before the first reconfinement shock. After this shock, the 
Mach number oscillates around the initial value and decreases, 
with slight increases in the expansion regions, until the head 
of the jet. Here, transonic and subsonic velocities result from 
the increase in temperature and decrease in velocity, as the flow 
crosses the reverse shock. The cocoon and the shell material are 
still in high overpressure with respect to the ambient by the end 
of the simulation, in a factor > 10 4 . The gas that forms the co- 
coon is initially fast and slightly supersonic, with velocities up to 
9 x 10 9 cms" 1 and Mach numbers up to 2, close to the head of the 
jet. Farther downstream, the backflow gets slower and subsonic, 
with velocities < 6 x 10 8 cms" 1 . 




1x10 19 2x10" 3x10 19 4x10 1s 
z (cm) 



Fig.A.l. Pressure along the jet axis as a function of distance 
from its injection point, at ~ 10 18 cm, as found in the hydrodi- 
namical simulations. Several conical shocks are present, due to 
the pressure balance with the surrounding cocoon: each time the 
jet pressure falls below that of the cocoon, a shock is formed, 
keeping the jet radius roughly constant until it reaches the re- 
verse shock. 
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Fig.A.2. Evolution of the self-similar parameter R = l\,lr\, as 
a function of time. After the pronounced initial increase, R re- 
mains between 2.5 and 2.7 for most of the simulated time. 
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Fig.A.3. Mass density (top) and pressure (bottom) maps resulting from hydrodynamical simulations. The simulations were per- 
formed using gjet = 3 x 10 35 erg s _1 , f M Q = 3x 10 4 yr and « ISM =0.3 cm -3 . 
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Fig. A.4. Lorentz factor (top) and Mach number (bottom) maps resulting from hydrodynamical simulations. The adopted parameters 
are the same as those of Fig. lA31 
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Fig. A.5. Pressure (left) and mass density (right) evolution in the shell (top) and cocoon (bottom) regions as a function of time. A 
fitting of the results is also shown for both variables (dotted line). This fit serves to estimate the simulation values at longer times. 
This extrapolation is strictly valid only if an homogeneous external medium and a constant injection energy rate are assumed. 



